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SUMMARY 


The effect of static tensile load on the thermal expansion of Gr/PI com- 
posite material was measured for seven different laminate configurations, A 
ccmiputer program was developed which implements laminate theory in a piecewise 
linear fashion to predict the coupled nonlinear thermcxnechanical behavior. 

Static tensile load significantly affected the thermal expansion charac- 
teristics of the laminates tested. This effect is attributed to a fiber- 
instability micromechanical behavior of the constituent materials. Analytical 
results correlated reasonably well with free thermal expansion tests (no load 
applied to the specimen) . However^ correlation was poor for tests with an 
applied load. 


INTRODUCTION 

Fiber-reinforced composite materials have been used extensively in aircraft 
secondary structure (refs. 1 and 2) during the past decade. Presently, many 
of the questions relating to the mechanical properties and fabrication tech- 
niques have been answered which has given the structural designer sufficient 
confidence for application to primary structure (refs. 3 to 5) . Future designs 
can take advantage of the inherent characteristics of composites to tailor the 
structural response to meet particular design requirements (ref. 6) . 

Before full advantage of the tailorability of composites can be achieved, 
more in-depth study must be made to understand the various thermomechanical 
behaviors. One such phenomenon is the effect of load on the thermal expansion. 

The effects of static tensile load on the coefficient of thermal expansion 
of mild steel and Invar were examined to a limited degree in reference 7. 
Resistance-type strain gages and a dead-load tensile test machine with a ther- 
mal chamber were used. A thermally calibrated material (Invar) was used to com- 
pensate for apparent thermal strain. In all tests the specimens were reported 
to reach a thermal and mechanical equilibrium prior to data being recorded. 

Using this test method, the authors reported a resolution on the order of 
10”^ (cm/cm) /K. 

A more recent study (ref. 8) was made of the effects of static compressive 
load on the thermal expansion of a [±45]g laminate composed of Celion 6000^ 
graphite fiber with PMR-1 5 polyimide resin. The results of that study showed 
that the static compressive load affected the thermal expansion of the composite 
material but, due to the limited number of tests and the reported quality of 
specimens tested, the magnitude of the effect could not be clearly established. 


^Celion 6000 is a registered trademark of the Celanese Corporation. 



For the study reported in reference 8 , resistance- type strain gages were 
used to measure strains and thermocouples were used to monitor the specimen 
tanperature. An IIT/Research Institute compression specimen was used in a 
thermal chamber with load applied by a displacement-controlled test machine. 
Compensation for apparent thermal strain was accomplished using quartz as a 
control material- All tests were performed using a continuous heating cycle 
at a constant heating rate of 83 K/hr. Data were recorded without insuring 
that the specimen had reached an equilibrium state. 

The present study examines the thermal expansion characteristics of Celion 
6000 graphite fiber with PMR-1 5 polyimide resinr between room temperature (RT) 
and 560 K. These tests were performed for free thermal expansion (no applied 
load) and applied static tensile load conditions. A computer program was 
developed which implements laminate theory in a piecewise linear fashion to pre- 
dict the coupled nonlinear thermanechanical iDehavior. Comparison of the analyti- 
cal and experimental results are made. Seven different laminates - [O 4 ] , [ 9 O 4 ] , 
[±30] gf [± 45 ] g, [ 0 / 90 ] gr and [ 90/-45/0/45] g - were tested. The effect of static 
tensile load on the thermal expansion characteristics of each laminate is com- 
pared with its free thermal expansion behavior. 

The use of trade names in this publication does not constitute an 
endorsement, either expressed or implied, by the National Aeronautics and 
Space Administration. 


SYMBOLS 


[A] 

CTE 

Ell r^22 
2 

{n} 

[Q]i 

RT 

T 

At 

ti 

,«22 


laminate stiffness matrix 
coefficient of thermal expansion 

tangent modulus in fiber and transverse-to-f iber direction in principal 
material coordinate system 

shear modulus in principal material coordinate system 
vector of normal stress resultants 

ith lamina constitutive matrix oriented in laminate reference frame 

room temperature 

temperature 

temperature increment 

thickness of ith lamina 

coefficients of thermal expansion in fiber and transverse-to-f iber 
direction in principal material coordinate system 
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ith lamina coefficient of thermal expansion oriented in laminate 
reference frame 

{a} vector of laminate coefficients of thermal expansion 

e strain 

{e} vector of midplane strains 

gAPP apparent thermal strain of gage 

thermal strain 
T 

e-ll thermal strain in fiber direction 

T 

£22 thermal strain transverse to fiber 

v-| 2 Poisson's ratio in principal material coordinate system 

a stress 

qULT ultimate stress 

T*| 2 shear stress in principal material coordinate system 


TEST SPECIMENS AND INSTRUMENTATION 

The experimental test matrix for this program is presented in table ^ . 
Seven different laminate configurations - [O4] , [9O4] , [+30] [±45]g, 

[±60] g, [ 0 / 90 ] g, and [ 90/-45/0/45] g - were tested with varying applied load 
conditions over a temperature range from RT to 560 K. All specimens were kept 
in a laboratory environment prior to testing. Based upon previous tests, mois- 
ture absorption was not considered a significant problem. The fiber volume 
fraction for all laminates tested was between 0.44 and 0.47. Nominal specimen 
dimensions are 20.3 cm long by 2.5 cm wide with a 12.7-cm gage length as 
depicted in figure 1 . Load introduction tabs were fabricated with Gr/PI 
identical to the test laminate to minimize any mismatch in thermal expansion. 
High-temperature strain gages (BLH FSM-1 2-35-SO) , placed axially back-to-back 
and transversely, were bonded to all specimens using a polyimide adhesive 
(BLH PLD-700) . The effects of gage transverse sensitivity were examined and 
found not to appreciably affect the measurement of the thermal expansion 
behavior of the laminates tested. 

The resistance -type foil strain gages used in this study are capable of 
operating in excess of 560 K. Electrical wires about 1.5 m long were silver 
soldered to the strain gages. For all tests, individual gages were connected 
to a full bridge using a three-wire system. 


A digital thermcxneter with type E thermocouple wire (Nickle chromium/ 
constantan) attached to the specimen was used to monitor the temperature. The 
temperature was displayed in increments of O.l^F to assist in determining when 
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thermal equilibrium was reached. Wire thermocouples were selected rather than 
other devices because of their greater versatility and temperature sensitivity 
and because they are not affected by coupled thermomechanical behavior of the 
test specimen. 

The calibrated specimen, as shown in figure 1 , used to compute the apparent 
thermal strain was made of titanium silicate. The specimen and detailed data 
on its expansion characteristics were supplied by the Corning Glass Works. The 
coefficient of thermal expansion (CTE) for titanium silicate is presented in 
figure 2. For all calculations in this study the mean value was used. Six 
strain gages, three back-to-back, were adhesively bonded to the titanium sili- 
cate specimen. Tests were conducted under the same thermal conditions as the 
composite specimens except for the application of load. The apparent thermal 
strain of the gage, as depicted in figure 3, was computed by taking the 

difference of the measured and calibrated responses of the titanium silicate. 

A titanium coupon specimen (fig. 1) 26.7 cm long and 2.5 cm wide was used 
to determine the applicability of the thermal expansion test method. Strain 
gages were oriented axially, transversely, and at 45° relative to the longitudi- 
nal axis of the specimen. Results fran this test are discussed later. 


TEST PROCEDURE 

The test procedure consists of the following sequential steps: 

(1) Align and attach test fixture to specimen. 

(2) Install fixture and specimen inside furnace to creep-frame load train 
and connect all instrumentation. 

(3) Apply static load. 

(4) Balance gages and record room- temperature strain readings. 

(5) Adjust set point temperature approximately 28 K above previous set 
point temperature. 

(6) Thermally condition specimen until equilibrium is reached 
(approximately 15 to 40 min) . 

(7) Record specimen temperature and strain. 

(8) If specimen temperature is below 560 K, repeat steps (5) to (7) . 

(9) Reduce temperature to approximately 450 K, allow specimen to reach 
equilibrium and record temperature and strain. (Step (9) is performed only for 
those tests where load was applied to the specimen for determining whether creep 
or plasticity occurred during the test.) 

The static tensile load was applied to the specimen through a mechanical 
advantage in a dead-load creep frame as depicted in figures 4 and 5. Specimen 


4 


grips and alignment apparatus, illustrated in figure 6, were used to introduce 
a uniform tensile load into the specimen. 

The magnitude of the maximum applied load for each different laminate was 
limited to the lesser of the load necessary to produce a strain of 0.004 at 
589 K or the load corresponding to the laminate initial nonlinear stress-strain 
response at 589 K. The 0.004 strain limit was employed because it corresponds 
to a typical composite-material design strain. 

A specimen was considered to be in equilibrium at a set point when the 
transient response of the temperature and strain dampened out and became in 
phase. At equilibrium the fluctuation in temperature was on the order of 
+0.2 K whereas the strain varied approximately ±3]ie. 

An examination of the effects of heating rate on the longitudinal expansion 
characteristics of a [ 9 O 4 ] Gr/PI composite laminate was performed. The effects 
of four heating profiles from RT to 560 K to RT, as depicted in figure 7, were 
examined. Three profiles had constant heating rates and the fourth a step pro- 
cedure. Based upon heat conduction laws, in the limit, the lower the heating 
rate, the nearer an equilibrium condition is achieved. The step heating pro- 
cedure (profile III) is capable of producing an equilibrium condition at each 

step, provided the soak time is sufficiently long. 

The thermal expansion of a [ 9 O 4 ] specimen was examined using each heating 
profile twice. The results shown in figure 8 correspond to only one test of 
each heating profile since good agreement between repetitive tests was obtained. 
The sharp contours of the heating profiles were not reproduced because of the 
thermal inertia of the test chamber and specimen. Upon examination of the ther- 
mal expansion results of the four heating profiles the following trends were 
obtained. A nonequilibrium condition existed for heating profiles I and II; 
this is evident from the reduced hysteresis between the results produced for 
heating profiles I, II, and IV. The step heating profile of III gave results 

similar to those of profile IV, depicting an equilibrium or near-equilibrium 

state. 


Based upon these results, a step heating profile was used on all subsequent 
tests. This heating profile, given sufficient soak time, will result in a state 
of equilibrium during data recording. Equilibrium was assumed to exist when the 
specimen temperature and the monitored strain readings became stabilized, gen- 
erally resulting in soak times between 15 and 40 min. 

A free thermal expansion test of a titanium alloy (Ti-6A1-4V) specimen was 
performed to establish the general validity of the test method. This test was 
performed between room temperature and 560 K. The titanium alloy was selected 
for its stable thermal expansion characteristics, the magnitude of the coeffi- 
cient of thermal expansion, and the availability of published data (ref. 9) . 

The thermal strain was measured and compared with published results as depicted 
in figure 9. Though scxne discrepancies exist between the measured and pub- 
lished results, less than 60 us, they are within nominal differences of mechan- 
ical properties between melts of the same material. This test was repeated and 
yielded consistent results. 
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ANALYSIS 


Classical laminate theory was used to ccmipare the predicted response with 
the measured thermomechanical response. A computer program was written to 
implement laminate theory in a piecewise linear manner to predict the overall 
nonlinear behavior. The equation governing the laminate response to mechanical 
load, neglecting bending, is 


{n} = [A(e,T)]{e} 


( 1 ) 


where 


M 

[a(€,T) 1 = ^ [Q(e,T)]iti (2) 

i=l 

The [q] ^ matrix is defined as the constitutive matrix of the ith lamina, 
relating to stress and strain, oriented in the laminate reference frame. The 
thermal expansion of the laminate is governed by 


= {a{a,e,T)> At 


(3) 


The vector of coefficients of thermal expansion for the laminate { 5 } is defined 
by 


M 

{a} = [A(e,T)]"^ 2! [Q(e,T]^ti 

i=l 


(4) 


where is the CTE of the ith lamina relative to the laminate reference 

frame. For a derivation of equation (4) the reader is referred to references 10 
and 1 1 . 

As a thermal or mechanical load is applied, in small increments, the lami- 
nate material properties are updated corresponding to the temperature and stress 
or strain state of each ply. Load increments on the order of 5 K and 1 percent 
of were found to provide converged solutions. 

The residual thermal stresses and any mechanically applied load are 
included in the analysis prior to the thermal expansion analysis. The residual 
thermal stresses are modeled by starting from the stress-free temperature (560 K) 
and reducing the temperature to room temperature. If a mechanical load is 
applied, then an analogous procedure is performed. An unsymmetric [0/90] 
laminate was heated to 589 K to determine the stress-free temperature. The 
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unsyiiiinetric laminate at room temperature is deformed due to the residual thermal 
stresses. When heated to 560 K the specimen became flatf thus indicating a 
stress-free state. 

The constitutive nonlinear material properties used in the analysis, E*| ] ^ 

T 

^22f and V-| 2 r are functions of strain and temperature whereas and 

T 

€22 are functions of stress and temperature. The constitutive properties were 
determined from the data presented in figures 10 to 13. 


CXDMPARISC»J OF THEORY AND EXPERIMENT 

The experimental results of the effect of static tensile load on thermal 
expansion are presented in figures 1 4 to 20. When repetitive tests were per- 
formed, bars denoting the nominal scatter band are used and the curve presented 
is the average of all tests. Each curve is a least-squares third-order poly- 
nominal representation of the experimental data. In each test nine data record- 
ings were made during the heating cycle at increments of approximately 28 K. 

The constitutive thermal expansion properties used in the analysis portion 

T T 

of the study, that is, e-|i and ^22^ developed using [O 4 ] , [ 9 O 4 ] , and 

[± 45 ] s laminates. The [±45]g laminate subjected to axial tension was used 

T T 

to determine the effect of shear stress on and ^ 22 * validity of 

these data is suspect because only a limited number of [±45]g specimens were 
available for testing and the resulting thermal expansion trends above 500 K 
were not expected. 

A static tensile stress produced a significant effect on the longitudinal 
thermal expansion of the [O 4 ] laminate as depicted in figure 14. For the free 
thermal expansion tests a small strain scatter band on the order of 0.0001 
resulted, whereas for the cases where load was applied the maximum scatter band 
was on the order of 0.0004. The thermal expansion curves presented in figure 14 
illustrate a characteristic "hump” between temperatures of 420 K and 500 K. The 
location of the hump shifts to lower temperatures with increasing applied load. 
This phenomenon may be attributed to the compressive residual thermal stresses 
in the fibers as a result of the differences in the coefficients of thermal 
expansion between fiber and matrix. 

Although explicit evidence cannot be provided, one hypothesis of the hump 
phenomenon is as follows: During the cool-down phase of the cure cycle, the 

mismatch in thermal expansion between fiber and matrix produces a compressive 
residual stress state in the fibers. These compressive stresses increase until 
a local fiber instability occurs which alters the apparent stiffness of the 
laminate and results in a hump in the thermal strain curve. This hypothesis 
is consistent with the results of fiber microbuckling studies reported in 
-^^rences 12 and 13. 
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The free thermal expansion test is the reverse process of the cool-down 
procedure. As the specimen is heated, the fibers tend to be straightened by 
the expansion of the matrix material. Once the fibers become straightened the 
matrix contributions are negligible and the overall expansion characteristics 
are representative of the fiber. For the case where a static load is applied 
to a specimen, the load tends to partially straighten the fibers. Upon heating 
the specimen, the hump occurs at a lower temperature because the fibers have 
been partially straightened by the applied load. 

As depicted in figure 15 the transverse thermal expansion characteristics 
were investigated using a unidirectional laminate with fibers oriented at 90® 
to the direction of the applied load. Load apparently had a lesser effect on 
the expansion characteristics than was the case in the [O 4 ] laminates. One 
important feature was noticed in that once load was applied to the specimen the 
magnitude of the load had little effect on the thermal expansion characteristics 
of the laminate. 

Based upon the results depicted in figure 1 6 the effects of shear stress 
T-| 2 oi' the longitudinal thermal expansion, in the fiber direction, was found 
to be negligible up to 500 K. At temperatures above 500 K plasticity effects 
contributed to the measured expansion of the laminate. Due to a limited number 
of specimens, repetitive tests were not performed; neither was the applied load 
reapplied at a lower level. The validity of the (±45] g laminate thermal expan- 
sion characteristics is suspected in light of the previously mentioned plastic- 
ity effects and results fron the (0/90]g and [90/-45/0/45] g laminates that are 
discussed herein. Further discussion of the effects of shear stress on a-ii is 
provided subsequently. 

As shown in figure 17, tensile load produced a significant effect on the 
thermal expansion of the [±30]g laminate. The experimental strain scatter was 
less than 0.0002, for specimens tested with and without applied load. The mag- 
nitude of the thermal expansion and the overall shape of the thermal strain 
curves for the [±30]g laminate are similar to those for the (O 4 ] laminate. The 
general shape of the thermal expansion curves was predictable although the 
agreement between experiment and analysis was inadequate. The measured thermal 
expansion of the laminate under load is greater than that of the unloaded lami- 
nate. The calculations depict the opposite trend, which could be a result of 
the effects of shear stress on otii and 0122 * In this analysis the effects 
of shear stress an and 0.22 have been neglected. 

As depicted in figure 1 8 an applied load produced a moderate effect, about 
20 percent, on the thermal expansion characteristics of a (±60]g laminate. At 
temperatures below 400 K the effect of an applied load was less than 2 percent, 
whereas above 400 K the difference increased to 20 percent at 500 K. The mea- 
sured free thermal expansion was greater than the expansion measured under 
applied load. This behavior was similar to that exhibited by the [ 9 O 4 ] lami- 
nate. Both of these laminates have matrix-dominant responses with minimal 
shear interaction. At temperatures below' 400 K good agreement, less than 
2 percent difference, between the experiment and analysis was obtained. At 
higher temperatures the measured free thermal expansion curve was not 
predictable. 
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Cross-ply [0/90]g and quasi- isotropic [90/-45/0/45] g laminates were 
tested with and without an applied load. Both laminates tested demonstrated 
that a static load had a significant effect on their thermal expansion behavior 
as illustrated in figures 19 and 20. The overall thermal expansion of these 
laminates is similar to that of other fiber-dominant laminates, such as the 
[ 04 ] and [±30]g laminates. The cross-ply and quasi-isotropic laminates strain- 
temperature curves have a characteristics hump similar to the [O 4 ] and [±30] g 
laminates although the hump occurs at a higher temperature for the cross-ply 
and quasi- isotropic laminates. This difference is suspected to be a result of 
the matrix-dominant response of the 90^ and 45^ plies. These matrix-dominant 
plies reduce the percent of 0 ^ plies which in effect reduces their contribution 
to the overall behavior of the laminate; hence, the characteristics of the hump 
are altered. Based upon the results of the [O 4 ], [0/90]g, and [90/-45/0/4S]g 
laminates, with decreasing percent of 0 ^ fibers the hump shifts to the right, 
that is, occurs at higher temperatures. 

The thermal expansion curve of the cross- ply laminate subjected to an 
applied load shifted below the free expansion curve analogous to the response 
of the [O 4 ] laminate. The quasi-isotropic laminate produced results similar 
to those of the [±30]g laminate. This vertical shift of the thermal- strain 
curve is attributed to the angle plies in the [±30]g and quasi- isotropic 
laminates. 

For both the cross- ply and quasi- isotropic laminates good agreement, dif- 
ferences of less than 10 percent, between experiment and analysis was obtained 
for the free thermal expansion tests. For the applied load conditions, agree- 
ment between e 3 q>eriment and analysis was poor. Furthermore, in the case of the 
quasi- isotropic laminate, the measured shift in the expansion curves due to load 
was the opposite of that predicted. 

Based upon laminate theory the free thermal expansion of [±45]g, [0/90]g, 
and [90/-45/0/45]g laminates should be identical. Close agreement in measured 
free thermal expansion was found for the [0/90]g and [90/-45/0/45] g laminates 
while the [±45]g did not agree although the l±4S]g and [0/90]g laminates were 
machined fron the same panel. The measured thermal expansion of the [±4S]g lami- 
nate is suspected to be in error and the shear stress effects on the longitudi- 
nal coefficient of thermal expansion a-]*! may be more significant than these 
measurements would indicate. 

If shear stress has a significant effect on 0^11 / then the inclusion of 
angle plies in fiber-dominant laminates could affect the thermal expansion 
trends. An attenpt was made to include the effects of shear stress t -|2 on 
the longitudinal thermal e 3 q>ansion ot-|i. Because the data for the [±45jg lami- 
nate were suspect, data relating shear stress and a-ii were generated such that 
when incorporated in the analysis, the thermal expansion characteristics of the 
[±30]g laminate subjected to load were predictable. The agreement between anal- 
ysis and experiment for the quasi- isotropic laminate under load was improved 
using these data even though the analysis predicted a decrease in thermal 
expansion under load rather than the increase measured. Further tests to deter- 
mine the effect of shear stress on a-ii and a more sophisticated analysis are 
probably required to adequately model the coupled behavior in laminates with 
angle plies. 
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OONCLUSI(»^S 


Based upon results obtained in this study the following conclusions can be 
drawn about the effects of applied load on the thermal expansion characteris- 
tics of Gr/PI composite laminates, 

1 • Static tensile load has an effect on the thermal expansion characteris- 
tics of the laminates tested. The magnitude of the effect is a function of the 
laminate orientation and, in general, the magnitude of the applied load. 

2. For the laminates tested, the [O 4 ] laminate was most strongly affected 

by the applied load. Fiber instability induced during cool down frcxn the stress- 
free condition near the cure temperature may be a significant factor. 

3. The test method used provided adequate resolution for depicting trends 
with a reasonable degree of scatter. 

4. Reasonable correlation was obtained between experiment and analysis, 
based upon classical laminate theory, for the free thermal expansion tests. 

Poor correlation was obtained for the applied load tests. A formulation taking 
into account the coupled material behavior in a more rigorous fashion may be 
beneficial. 

5. Shear stress Ti 2 could have a significant effect on the longitudinal 
thermal expansion coefficient a*|i although the magnitude of its effect was 
not conclusively established. Further experimental and analytical studies 
are required to account for shear stress effects. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 1 , 1 981 
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TABI£ 1 EXPERIMENTAL TEST MATRIX 


Specimen 

Laminate 

layup 

Applied 
stress, MPa 

Specimen 

Laminate 

layup 

1 

[O4] 

0 

1 4 

[904] 

2 

[O4] 

0 

1 5 

[±60]g 

3 

[O4I 

287.0 

1 6 

[±60]g 

4 

10,1 

287.0 

1 7 

[±60]g 

5 

[O4] 

403.0 

18 

[0/90] g 

6 

[O4] 

403.0 

1 9 

[0/90lg 

7 

[9O4] 

0 

20 

[0/90] g 

8 

[9O4] 

0 

21 

[±30]g 

9 

[9O4] 

8.6 

22 

[±30]g 

10 

[9O4] 

8.6 

23 

[90/-45/0/45]g 

11 

[9O4] 

1 3.4 

24 

[90/-45/0/45]g 

1 2 

[9O4] 

13.4 

25 

[90/-45/0/45]g 

13 

[9O4] 

18.2 

26 

[± 45 ] g 


Applied 
stress, MPa 

18.2 

0 

0, 16.2 
0 

0, 294 
0, 294 
0, 294 
0, 105 
0, 105 
0, 158 
0 

1 58 


0, 9.9 






500 


550 


600 


TEMPERATURE, K 


Figure 2.- Coefficient of thermal expansion of titanium silicate. 























Figure 7.- Heating profiles investigated to determine effects on apparent strain 
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Figure 10.- [ 04 ] temperature- dependent stress-strain curves of Gr/PI 

for determining E-]i(a,T). 
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300 K 


Figure 11 *- [904] temperature- dependent stress-strain curves of Gr/PI 

for determining E22 (^fT) • 
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SHEAR STRAIN 

Figure 13.- Temperature-dependent in-plane shear stress-strain response of 
r±45]g laminate for determining G*| 2 ('^fT). 
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TEMPERATURE, K 

Figure 14.- Measured longitudinal thermal strain of [O 4 ] Gr/PI laminate 
as function of temperature and static tensile stress. 
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5 .- Measured longitudinal thermal strain of [904] Gr/PI laminate 
as function of temperature and static tensile stress. 
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Figure 16.- Thermal strain in fiber direction in [±45]g laminate as function 
of tanperature and applied static shear stress Tio. 



Figure 17.- Measured and calculated longitudinal thermal strain of f±30] 
Gr/PI laminate as function of tenperature and static tensile stress. 
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Figure 18.- Measured and calculated longitudinal thermal strain of r±60]g 
Gr/PI laminate as function of temperature and static tensile stress. 
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Figure 19.- Measured and calculated longitudinal thermal strain of [0/90] 
Gr/Pl laminate as function of tenperature and static tensile stress. 
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